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Abstract. - The entanglement-sharing properties of an infinite spin-chain are studied when 
CO ' the state of the chain is a pure, translation-invariant state with a matrix- product structure [1]. 

^ . We study the entanglement properties of such states by means of their finitely correlated struc- 

ture [2]. These states are recursively constructed by means of an auxiliary density matrix p on 
a matrix algebra B and a completely positive map E : ^ ® S — > B, where A is the spin 2x2 
matrix algebra. General structural results for the infinite chain are therefore obtained by ex- 
plicit calculations in (finite) matrix algebras. In particular, we study not only the entanglement 
, shared by nearest-neighbours, but also, differently from previous works [3], the entanglement 

■ shared between connected regions of the spin-chain. This range of possible applications is il- 

lustrated and the maximal concurrence C = [4] for the entanglement of connected regions 
can actually be reached. 
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^ . The recent developments in quantum information and computation theory are witnessing 
an increasing interest in the entanglement properties of multi-qubit systems Hke quantum 
^ ,. spin-chains, correlated electrons and interacting bosons. These systems had so far been stud- 
' ied in condensed matter physics especially in relations to their critical behaviour in phase- 
. transitions [5]. With its emergence as a precious resource for quantum computation and 
' communication tasks, entanglement is now investigated in condensed matter physics, too; on 
one hand to detect, extract and manipulate it, on the other hand to understand its role in 
solid state phenomena [6-13]. 

The entanglement properties of the class of states of quantum spin chains known as Matrix- 
Product States (MPS) [1], have been proved powerful tools in density matrix renormalization 
groups techniques [14] and in universal quantum computation to show the equivalence of 
teleportation-based and one-way quantum computation [15]. They can be also achieved by 
sequential generation [16]. Actually, these states were originally introduced in [2] as Finitely 
Correlated States (FCS) and their properties were there studied in great detail; in particular, it 
was showed that FCS can be implemented as ground states of Hamiltonians with short range 
interaction (see also [15]). This motivates the study of entanglement in FCS both from the 
point of view of the large variety of scaling laws provided and from the point of view of the 
possible physical applications. 
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Two main different approaches tiave emerged: one deaiing witti tlie scaiing beliaviour of 
entangiement-slraring between various subsets of localized spins [17, 18, 20, 21], the other one 
studying how much entanglement can be localized on two distant spins by measuring the 
others [19,22]. In this Letter we choose the first approach and study the entanglement of FCS 
by means of their recursive structure [2], that is by studying the specific completely positive 
maps between finite dimensional algebras and finite dimensional density matrices on which 
it is based. This provides a rich phenomenology of total translation invariant states over 
infinite spin chains whose entanglement properties are nevertheless determined by, and can 
consequently be studied in a finite-dimensional setting. Thus, we are able to compare nearest 
neighbours with non-nearest neighbours entanglement or, more generally, to analytically study 
the entanglement between different subsets of lattice points. In particular, we shall give a 
general necessary condition for one spin being entangled with a subset of others and prove 
that in some cases this condition is also sufficient. Finally, we show that entanglement sharing 
can achieve its maximum in the sense of [4]. 

Construction of translation invariant FCS: We will denote by Az an infinite spin- 
chain, the spins at sites i G Z being described by the algebra {A)i — M2 of 2 x 2 complex 
matrices. The infinite algebra Az arises as a suitable limit of the local tensor-product algebras 
A[-n,n\ •= ®j=_n(-^)j- ^'^y State CO Over Az is specified by density matrices P[i.„] defining 
the action of w as an expectation over local operators A^i ^] G -4[i „]: 

w(A[i,„]) = Tr[i^„](^P[i^„]A[i^„]) . (1) 

The /9[i^„]'s must satisfy the compatibility conditions 

Tr„+i(p[i,„+i]A[i_„] (g) l„+i^= Tr(^P[i_„]yl[i_„]) , (2) 
whereas, translation-invariance requires 

Tr„+i(p[i,„+i]li «) A[2,„+i])= Tr(^P[i^„]A[i^„]^ . (3) 

The class of translation-invariant FCS over Az is defined by a triple {B, p, E) where B is a. bxb 
matrix algebra B, p ^ B a, density matrix and : A® B ^ B a. completely positive unital 
map, which in Kraus-Stinespring form reads 

¥.{A(DB)^^Vj{A(dB)v] , V, : £^ d) C!' ^ , (4) 
i 

with A G A and B G B. Unitality means that identities are preserved: E(l^ (X) Ig) = Ig. 

Let ¥,^^'>{A) :— E(A (g) Ig); this defines a completely positive map from A into B. Analo- 
gously, the recursive compositions E^") := E o ^id^ (g E^^""^^^ are completely positive maps 
from .4[i „] into B. Setting 

Tr(p[i,„] A[i,„]) :=TrB(pE(")(A[i,„])) , (5) 

the r.h.s. recursively defines local density matrices P[i,n] over and a total state to on 

Az [2]. Further, Trg (^pE(^l^(g)B^ )^ = TrB(pB), VB G S yields translation-invariance [2]. 

Concretely, we choose B = M2 and E in Q with just one Kraus operator V : C^0C^ i~> C^. 
This is such that V\(f>^ (g V) = v^\'^p), V^^/j) = and 



E 



([(/),) (0, 1 (g , BeB, (6) 
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where \4>i.2) G are orthonormal and vi^2 are 2x2 matrices satisfying 

2 

viv\ + V2vl — 1b (unitality) , wj p Vj = p (translation invariance) . (7) 

If there exists a unique p fulfilling the previous condition, the resulting translation-invariant 
FCS are pure states over Az [2], namely they cannot be decomposed as mixtures of other 
states. These pure states can be interpreted as ground states for appropriate constructed 
Hamiltonians of finite range, but higher correlations [2] . 
Consider Q with n = 2 and ^[1.2] = Ai (E) A2, then 

Tr(/9[i.2]Ai® A2) :=Tre(pE(Ai(gE(A2(8)lB)) . (8) 

Using the properties of the trace-operation, the action of E becomes the action of its dual 
map F onto the state p e B: Trjsip E{A ® B)) = Tr^,ge(F(/9) A (K) B). This provides a state 
PAm ■■= F(p) = V^pVoiiA(g>B: 



V^pVi V^pV2 



The r.h.s. of ||SJ) reads Tr_4,g,e (^PA(g>B Ai ® ^E^^2 ® Ig^ ) ) ' ^■^ turning id^ ® E into its dual, 
nearest-neighbours states arise as pi2 '■— P[i.2] = Trglid^ (g) F(p^,g,B) j and read 



2 

ijlm—1 


/ Rlijl 

\R2ijl 




-R1112 


-R1121 


-R2111 


-R2112 


-R2121 


-R1211 


-R1212 


^1221 


\-R2211 


-R2212 


^2221 



Rlij2 
R2ij2 

^1122\ 

^2122 

i?1222 



(10) 



where Rijim = Tr(?;|z)j p viv„i), while general local density matrices are given by 

Pll,n] PVt) , (11) 

s,t 

where \4)s) = \4>si (t>s2 ® ■ ■ ■ 4>s„), vt := wti • • • 

We shall now relate the entanglement of pa(^b to that of the spin at site 1 with those 
at sites in [p,n\, p > I. This means investigating the entanglement of the restricted state 
CO \{A)i A[p,n]- By the previous construction, the restriction amounts to the expectations 



uj(^Ai 



l[2.p-l] ®Ap(g,--- An^ = TTA(S)B[pA(g,B ^1 ® <G{Ap (g, ■ ■ ■ A„)) , (12) 



where l[2,p-i] — (1)2 (8) • • • (l)p-i and G is a completely positive map from A[p^n] into B. 
Notice that while pi2 is a nearest-neighbours state, PA(8)B encodes the entanglement of the 
spin at site 1 with any subset of spins -4[p „], p > I, after being embedded into B = A = 
M2 by G. The advantage of the abstract structure presented above emerges in that the 
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total state over the chain is determined by the triple (S,E, p), so that properties like the 
maximal entanglement shared by nearest neighbours in a translation invariant chain, and, 
more generally, the entanglement between different subsets of lattice points, can be studied 
by means of that triple only. 

Distribution of the entanglement along the chain: As a measure of the entanglement 
of a state 1/12 over the tensor product algebra Afi ^2 , we shall use the entanglement of 
formation [25] 

Em.^M. (^12) = inf . {Y, \Afi) I , (13) 

where 1^12 = ^j^vi ^ convex decomposition into pure states and s{^\^ \^\^ the von 
Neumann entropy of their restrictions to the subalgebra M\ , namely of the state obtained by 
tracing over the second factor, the density matrix representing the expectation functional V\2- 
It turns out that 

E{A)x®{A)p^'^) ^ -^(^)i«>^[p,„]('^) ^ Ea®b{.Pa®b) ■ (14) 

The proof of this fact follows from a slight generalization of an argument in [24]. Let As 
be another algebra, v\-i a state over M\ ® Afs and T : N'2 ^ A/a a unital (Tllf^f.^) = 1^.^) 
completely positive map. Then, V12 = vi^ o (idj^^ ®T) is a state over A/i ® Af2, where o 
denotes the composition of the expectation functional 1^13 with the completely positive unital 
map id^Aj (g) F. As the decompositions of 1^12 induced by those of vis are not all possible ones, 
it follows that £'A/'i(giA/2 (^^12) is bounded from above by 

inf . f ^ A, S (v{^ o (idAA, ® r) [AAi ) } . (15) 

From the unitality of F, the states o {idj^^ F) restricted to JVi, that is evaluating mean 
values of operators of the form A^i (g) I7V2 , coincide with the restrictions of i^l^ themselves. 
Thus, E^^0^^{i/i2) < -EaAi«iAA3(j^13)- The first inequahty in the hierarchy ifT^ follows by 
taking A/i = N2 — (-4)p, = -A^[p,n] and as F the natural embedding of [A)p into 

A[p^n]- The second one, by choosing A/*! = (^)i, J^2 = -4[p.„], Nz — B and F G : A[p^n] ^ B; 
in fact, u) \{A)i ® -4[p,n] — PA®B ° id^ ® G, see also Eq. ((T^ . Thus, the general necessary 
condition follows that, for (^)i to be entangled with ^[p,„], p > 1, the state PA®B rnust 
be entangled over A® B. Translation-invariance makes this necessary condition independent 
of rigid shifts of the two algebras; also, though the next examples go in the right direction, 
unfortunately, there is so far no general argument for its sufficiency namely that entanglement 
between A and B should imply entanglement between A and some yl[p.„] . 

Examples: Because A^ B = M2, we study the concurrence [26] of which the entanglement 
of formation is a monotonically increasing function. The concurrence of p is given by C{p) — 
max{0, Ai — A2 — A3 — A4}, where Aj are the square roots of the eigenvalues in decreasing order 
of the matrix pp where p — [cry ® cry)p* {o'y ® cry) and p* denotes complex conjugation in the 
standard basis. In terms of concurrence, inequality H14|l reads 

c(uj\{A)i®{A)p) <c(Lj\{A)i®{A)^^r.^ <c(^PAm) ■ (16) 

As a first finitely correlated structure, let us choose vi ~ o) ' ^ (^s ' ^^^^ 
c — cos 95 and s = sin 1^9. Then the conditions Q are satisfied by /? = 1/2 f ^ ^ j and. 
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Fig. 1 - (a) Ca^b, green; C12, red; abscissa a = siri(2iy3). 



(b) Ca^b, green; C12, red; abscissa (p. 



using the state of three adjacent qubits exphcitly reads 
9 X! ® ^^ij ' where 

17n 



■2 2c's 

^2 



C 



2^2 



whereas the state pa^b amounts to 



s 



2c\s^ 
^2 



21 



cs 2c^s 
2cs^ cs 



I 



PA®B 



CS 
2c2s2 



\2c 



cs 
s" 
2cs^ 

2s2c2 



2c2s2 

2cs3 
cs 



2c3s\ 

2c2s2 



cs 

^2 



From (|17|l . tracing over site 2 and 3, one gets the one-site qubit-states pi 
while tracing over site 3 we get the following two-sites states, 



P12 



1 

4c2s2 





( 


2c2s2 


2c2s2 


4c4s2\ 


1 


2c2s2 


s2 


4c2s4 


2c2s2 


2 


2c2s2 


4c2s4 


s2 


2c2s2 




i^4c4s2 


2s2c2 


2c2s2 





(17) 



(18) 



(19) 



4c2s2 
1 



(20) 



and tracing over site 2 gives P13 = ^^^n (8) -|- 2^22 "22- 

The latter state is clearly separable, there is no entanglement between site 1 and site 3, 
i.e. no next nearest neighbour entanglement. While partially transposing pi2, i.e. Ac'^s^ ^ 
Ad^s^ the positivity is surely lost for < 1^9 < 7r/4, thus we have entanglement between 
nearest neighbours. The same happens to the state of one site with the rest pA®B, namely, 
if s^ — 4c^ < 0, that is for < tp < 7r/4, this state is surely not positive under partial 
transposition. 

A finer picture can be obtained by looking at the concurrences: Ca®b '■= C{pA(giB) = 
I sin2(/?cos(2(y9)| vanishes sA ip — 7r/4, 57r/4, while C12 ■— C{pi2), though explicitly computable, 
is not as easily readable. Their behaviours as functions of a = sin(2i^) are reported in Figure 
1(a) which shows firstly that, in agreement with (|15|l . Ca^b ^ ^12 and secondly that pa0B is 
entangled whenever pi2 is entangled, namely for all sm{2ip) 7^ 0,±1. 

Another finitely correlated state is given by choosing 



(21) 
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with s = sinip and c = cos 1^9. Conditions give 



y 

1-x 



+ 2s2 



(a— l)sc 



(22) 



With this choice (^2)^ = and all terms in (|10|l with two adjacent V2 or vanish; therefore, 
independently of a, s in (|22|l . nearest-neighbours states are of the form 



P12 



/I - 27 a a 0\ 

a 7^0 

a /3 7 

V 000/ 



(23) 



and have concurrence C12 = 2|/3|. The entanglement of formation is maximal when 

2(1 -a2)s2c2 



C12 = 2Tr(u|u|pw2Wi) 



(l-a)2c2 + 2s2 



(24) 



is maximal, that is when a = \/2 — 1, cos(2(/)i) = \/2 — 1 and /? = ^(v^ — 1), whence 
C12 = \/2 - 1 = 0.41. On the other hand, the concurrence of PA®B amounts to 



(l-a)2c2 + 2s2 



(25) 



In Figure 1(b), a is set V2 " 1 and I|24l - I25|) are plotted against (p £ [0, 27r]: again agreement 
is shown with the general monotonicity of concurrence expressed by (I16|l . Namely, nearest- 
neighbours are entangled if the state of one site with the rest PA^B is entangled, but we do not 
know if entanglement of one site with rest always implies nearest neighbours entanglement. 
We believe this to be peculiar of our choice of B: since it is two dimensional both pi2 and 
PA<»B have rank 2 and separable density matrices of rank 2 have measure. 

When three qubits share equal entanglement, then the maximum entanglement of one 
qubit with each of the two other qubits is the sum of the squares of the two concurrences 
which is less than or equal to one [4]. Thus nearest neighbour concurrence cannot exceed 
I/V2, but it is an open problem whether this upper bound can be achieved by an entangled 
chain. In [3] an entangled chain is considered in a translational-invariant state with nearest- 
neighbours states as in l|23|) but with a = and an upper bound for shared entanglement 
corresponding to a greater concurrence C12 = 0.434467. 

The translation invariant state over the system A[-oo,-i] '8)^o®-^[i,oo] can be defined with 
(5 (g) 1 ig) G mapping the system into B ^ A ^ B. Now we have a effectively three qubit state 
sharing equal amount of entanglement and due to the considerations in [4] the upper limit for 
the concurrence is given by In the following example, we show how the maximum can 

indeed be reached by a finitely correlated structure. We set vi as in l|21|) . V2 = \/T~~a2 ^ ^ 



and 



a 



1 
2 \ a2 



then pAm 



ac 






2 

a 

as 




ac 



^/T^ 













,(26) 
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whence C(pa^b) = o — simp which attains its maximum l/\/2 at (p = tt/2 and a = 

1 + 

l/\/3- This state differs from the one in Eq. (|23|l . Therefore optimizing the entanglement of 
Ao (8) ^[i.oo] reduces the entanglement of Aq^:^ Ai- 

Conclusions: We studied the entanglement properties of FCS over infinite quantum spin 
chains by means of their recursive structure and illustrated some possible behaviors by exam- 
ples with 2x2 matrices. Investigation of higher dimensional contexts as the AKLT model [23] 
and comparison with different approaches [17] is in progress. 

Acknowledgment: B.C. Hiesmayr acknowledges EURIDICE HPRN-CT-2002-00311. 
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